Spin-Hall conductivity due to Rashba spin-orbit interaction in disordered systems 
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We consider the spin-Hall current in a disordered two-dimensional electron gas in the presence of 
Rashba spin-orbit interaction. We derive a generalized Kubo-Greenwood formula for the spin-Hall 
conductivity and evaluate it in an systematic way using standard diagrammatic techniques for 
disordered systems. We find that in the diffusive regime both Boltzmann and the weak localization 
contributions to are of the same order and vanish in the zero frequency limit. We show that 
the uniform spin current is given by the total time derivative of the magnetization from which 
we can conclude that the spin current vanishes exactly in the stationary limit. This conclusion is 
valid for arbitrary spin-independent disorder, external electric field strength, and also for interacting 
electrons. 
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I. INTRODUCTION 



The idea of using the spin degrees of freedom of elec- 
trons in semiconductor systems instead of their charge 
has attracted wide interest in recent yearsiiS*^. Of par- 
ticular interest are systems with spin-orbit interaction 
as this allows access to spin via charge and vice versa. 
A number of effects in such systems have been studied 

thcOrcticall3T5i6i7i8i9il0,llil2il3,14,15,16,17,18,19,20,21,22^ 

experimentalljiiiSiSiiiSi^, and most recently by Kato et 
alm^ who report the first observation of spin Hall effect 
in GaAs and InGaAs samples of the extrinsic type^. In 
particular, there has been considerable interest in Rashba 
spin-orbit interaction^^ in two-dimensional electron sys- 
tems since this type of spin-charge coupling can be con- 
trolled by electrical gates. Moreover, it has been pointed 
out by Sinova et alM' that clean (without impurities) 
systems with spin-orbit interaction carry spin currents 
and show a spin-Hall effect: applying an electric field 
Ex in x-direction generates a spin current in perpendic- 
ular y-direction. In a number of subsequent papers the 
corresponding spin-Hall conductivity cr,y^ has been cal- 
culated, and there seems to be agreement now that the 
spin-Hall conductivity at zero frequency vanishes in an 
infinite systemiSiiSiSiiSLSSiSS. In this paper we present a 
systematic calculation of tr^^, by using well-known per- 
turbative techniques for disordered systems^. The sys- 
tematic expansion is performed in powers of the small 
parameter 1/ppl, with pp the Fermi momentum and I 
the mean free path. In addition to confirming the re- 
sults of earlier worlsiSiiSiSLSLSS, which were obtained in 
the Boltzmann (semiclassical) regime, we go beyond this 
limit and also include the weak localization correction 
diagrams. We show that this latter contribution is of 
the same order in l/ppl and Rashba amplitude as the 
Boltzmann value and vanishes as well at zero frequency. 
This demonstration requires us to consider the various 
contributions to the Hikami box in the weak localization 
contribution as well as the rcnormalization of the current 
vertices. 



We then go on and show that the vanishing of the spin 
current is expected on rather general grounds and under 
general conditions, such as arbitrary (spin-independent) 
disorder, electric field strength, and for interacting elec- 
trons. Crucial for our argument are two observations. 
First, the uniform spin current is given by the total time 
derivative of the magnetization^, which, in turn, pro- 
vides a straightforward interpretation of the spin current 
and its observable effect in terms of magnetization and its 
change in time. Second, a system driven by an external 
force which is constant in time normally reaches a steady 
state, implying then that the spin current vanishes under 
rather general conditions. 



II. MODEL SYSTEM 

We consider non-interacting electrons of mass m and 
charge e moving in a disordered two-dimensional system 
in the presence of Rashba spin-orbit interaction with am- 
plitude c^iiiSiS^. This system is described by the Hamil- 
tonian (we set h = 1) 



H=^+a{&'py-a'^Px) + U{r), (1) 

where p denotes the momentum operator, and are 
the Pauli matrices. Further, U (r) is a short-ranged 
disorder potential with the property that C/(r)C/(r') = 
(mT)~^(5(r— r '), where the overbar indicates average over 
the disorder configuration, and r the mean free time be- 
tween collisions. 

The spin-orbit term in the Hamiltonian |^ changes 
the expression for the charge-current density operator (in 
standard second quantization notation)^: 



2m 



—i'^ ( A - 
mc 



(2) 

where together with an electromagnetic vector potential 
A (coupling to the charge e) a fictitious spin-dependent 



2 



vector potential A is introdueed: 

A = — (Q!mc/e)(— (T^, (T"'^, 0), 



(3) 



where c is the speed of hght. On the other hand, 
the spin-current density operator (z-component of spin) 
remains the same as in the case without spin-orbit 
interactionSSiSi: 



j^-(r) = - ^ TA^fT^W^ - (VV3^)<T^?/' 



2ie 



AV;^(r)(7^?/>(r) 



(4) 



III. SPIN-HALL CURRENT 



Considering now a hnear response regime for the spin 
current in the presence of an electric field E^, we derive 
a generalized Kubo-Greenwood formula for the spin-Hall 
conductivity. In the Keldysh approachS^ the spin current 
can be expressed as 



(5) 



2m 



°° dE 

Tr 

27r 



^a) ^Gk{E,E-u;)}, 



where Gk is the Keldysh component of the 2x2 matrix 
Green function; Tr stands for the trace in both momen- 
tum and spin space. Due to the Pauli matrix the 
diamagnetic term cx A gives a vanishing contribution^ 
to the spin current in cq. ISJ . Then in first order pertur- 
bation in A we find 



SGK{E,E-Lo)^{hE-hE-^) 



e 
mc 



A) A(w) 



G 



E-u 



(6) 



where He = tanh{E/2kBT), with fceT < Ep the tem- 
perature and E-p the Fermi energy, and where G^y^ = 

{E — H + Ep ± iO+)~^ are the retarded and advanced 
Green functions, respectively. From eqs. (O and © we 
obtain a generalized Kubo-Greenwood formula for the 
spin-Hall conductivity^. 



Tr 



-AJG 



(7) 



where E,u; E-p. In deriving lO we have assumed that 
the electric field Ex{u>) = iLuAx{Lu)/c is applied along the 
x-dircction producing a perpendicular spin current along 
j/-direction, i.e.. 



(8) 



Next, after averaging over the random disorder, the 
spin-dependent Green functions Gr/a become diag- 
onal in momentum representation^, (p|Gp^^|p') = 



(2^)2Gf,A(p)^(p-p'), with 



G 



R/A 



sM 



2 ^^^E-ap)~sap±^' 



(9) 



where M = {pya^ - Px^^)/p with M'^ = a° = 1, £_{p) = 
p^ /2m — Ep. The expression Q has been obtained in 
the self-consistent Born approximation. 



G 




G a 




FIG. 1: Diagrams for the spin-Hall conductivity in the zero- 
loop approximation: The right diagram represents the vertex 
correction. Dashed lines denote averaging over the impurities. 

The diagrams of the standard perturbative approaches 
can be generated from Q by expanding in the disorder 
potential U and "dressing" the Green function lines with 
diffusons and Cooperons. In this formalism22i2&, the re- 
sulting diagrams can be estimated according to the num- 
ber of loops, composed by Cooperon and diffuson lines 
(loop expansion). In the following sections we will ex- 
pand CT^j. in powers of l/ppl, where / = vpT denotes the 



mean free path, and evaluate it up to order 



that 



is, we will neglect all terms with n > 1 in the expansion 



IE 



(ppiy 



ppi > 1. 



(10) 



Thus in our diagrammatic loop expansion we have to con- 
sider two classes of diagrams: (i) diagrams with no loops 
(semiclassical approximation), and (ii) diagrams with one 
loop (weak localization corrections). This is so because 
of the non-standard situation arising from the spin de- 
pendence of the charge current vertex in lO, which can 
be written as 



e ~ 

Px = ppn 

c 



X ~1~ {^x 



PFflx 



(11) 



where pp is the Fermi momentum, and n = p/p is an 
operator of the direction of momentum p. The contri- 
bution in the brackets on the right-hand side of is a 
correction of order {ppl)~^ ^ 1 to the main term ppfix- 
Thus, in zero-loop approximation (sec Sec. lIVp this cor- 
rection leads to a contribution which is of the same order 
as the one coming from ppfix in the one-loop approxi- 
mation (i.e. weak localization correction - see Sec. EJ- 
Hence, both of them are expected to give contributions 
of the same order in 1 /ppl to the spin-Hall conductivity. 
The spin current vertex UzPy, on the other hand, can al- 
ways be approximated by dzPpfiy in the order considered 
in this work. 
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IV. ZERO-LOOP APPROXIMATION 

In zero-loop approximation we need to retain only 
the diagrams depicted in Fig. ^ with no loops. This 
approximation corresponds to the semiclassical (Boltz- 
mann) limit. To calculate the first diagram, we simply 
substitute Ga/r hi (0) by the averaged Green functions 
One can see that the main term of the charge cur- 
rent vertex does not contribute in the first diagram. 
The remaining terms in the brackets in give 



A = 1 - 



ILUT, 



(12) 



where x = 2apFT is the dimensionless spin-orbit param- 
eter which measures the spin rotation between two con- 
secutive impurity collisions. For w = this expression 
agrees with the qualitative result of Ref>i^ in the limit of 
small spin-orbit interaction ma^r <C 1, as well as in the 
ballistic limit t oo. 

Next we turn to the vertex correction, described by 
the second diagram in Fig. ^ Again, it is crucial^ to 
retain the entire correction of the charge current vertex 
(tTT|l (which amounts to include curvature effects of the 
Fermi surface), since it will turn out now that this cor- 
rection term leads to an exact cancellation of CTyi°^ given 
in eq. ifT^ . We note that the scattering off impurities is 
spin-independent, so that dashed lines do not alter the 
spin direction which simplifies the calculation consider- 
ably. For further evaluation it is convenient to apply the 
following identity for the Pauli matrices^ 



'0 ~Q 

SlS2 S3S4 



1 ^ 



S4S2 



(13) 



to every dashed line in the vertex diagram Fig. ^ with 
a, /3 = 0, . . . , 3, and where, again, is the unity matrix. 
This trick reduces the vertex diagram to a geometric se- 
ries consisting of terms where the momentum and spin 
summation can be performed separately. These terms 
are given by 



TrKGf^(p)a^Gf-(p)], (14) 



2mT 



, 2 

where from now on Tr stands for / Tropin. The off- 
diagonal elements of X'^^ vanish, while the diagonal ones 
are given by 



~ A' "^^+A^^ 



-1^11 _ v22 _ 1 
— — — 



1 



A a;2+A2 



(15) 



by two "bubbles" with a diffuson wavy line in between: 




27rm2 



Py r<E 



G-(p)a^Gf--(p) 



,^MM Tr [a^Gf (p) (p, - ^i,) Gf-'^(p)] , (16) 



where D^^ denotes the diffuson 

1 1 



2mT 1 



(17) 



with X!^^ given in \ib\ . Only the term with /i = 2 con- 
tributes to (|16|l . The two bubbles in \i>6\ are explicitly 
given by 



Tr 



— — GR(p)cr 



^(P) 



xmrpF /2 
.t2 + A2 ' 



(18) 



and 



Tr 



a^Gf (p) (p. 



Gf--(p) 



A a;2 + A2 



(19) 



Inserting these expressions into H16(l , we finally obtain for 
the vertex correction of the spin-Hall conductivity 



x2/2 



in 



A2 {l-2iu;T) 2^ -iurX^' 



(20) 



Adding now eqs. (|12|l and (|20|) we obtain the spin-Hall 
conductivity in zero-loop approximation. 



yx\'~) "yx / ' "yx 

2 ^ a;2/2 



Stt a;2 + A2 



(1 - 2icjr)a;2/2 - itorX^ 



(21) 



which is the same result as obtained in RefiiS (see 
eq. (25) therein) by the equation of motion approach. 
Thus, in zero-loop approximation the two diagrams in 
Fig. n] cancel each other, so that the leading contribu- 
tion to the spin-Hall conductivity vanishes in the limit 
= 0, in agreement withi2ii2i2ii2L22i2S. A similar can- 
celation has been observed also for the case of the elec- 
tric conductivitji^. In both these cases the "anomalous" 
term cx A in the charge current operator Q is cancelled 
by the vertex correction at w = 0. Thus, the second di- 
agram in Fig. ^ can be interpreted as a renormalization 
of the charge current vertex which results in the cancel- 
lation of the "anomalous" term oc A^- 



V. WEAK LOCALIZATION CORRECTION 



The expression for the vertex diagram can be split into 
three parts, and the vertex diagram can be represented 



At the end of Sec. lIIII we argued that in order to reach 
a given accuracy, we have to include the contribution 
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For q = 0, the three contributions to the Hikami box 
in eq. (|25|l become 

Vr = Tr {L(p)a^ [a^i?(-p)]^} , (27) 



FIG. 2: Weak localization diagrams. The wavy line represents 
the Cooperon 12311 . Each diagram contains one spin current 
vertex (left) and one charge current vertex (right). 



of one-loop diagrams in our calculation of ct^^.. In this 
section we demonstrate this now by calculating the weak 
localization correction to ct^^ explicitly. 

Similar to the case of the charge conductivitji^, the 
weak localization correction to the spin-Hall conductivity 
is given by three diagrams depicted in Fig. |21 For sim- 
plicity, we calculate them at zero frequency uj = only; 
to simplify notations we omit the energy superscript of 
Green functions, assuming G^jx = G^^^ below. 

Using H13|) . we may separate the spin indices of a 
Hikami box from those of the Cooperon, like we did above 
in (|16|) . Then the total contribution of the three diagrams 
in Fig. 12 can be written as 



27rm2 ^ J (27r)2 



E 



C77 (q)v-77 (q)^ (22) 



where C'^'' (q) are the matrix elements of the Cooperon 
defined as 



C^"^ (q) 



1 



2mT 



^c(q) 
l-Xc(q) 



,7,7' = 0...3, (23) 



77' 



2mT 



TrKGR(p)a'^(GA(q-p))^]. (24) 



In eq. H22() . V^'^ (q) denotes the Hikami box: 
y77'(q) ^ y77'(q) + ^.77'(q) + ^77' (q). 



(25) 



Each of the three terms in eq. H25() corresponds to one di- 
agram in Fig.[21plus the corresponding diagram in Fig.|31 
(see eq. below). 

In the rest of this section we assume a:;^ ^ 1 This 
permits us to approximate H22() by 

27rm2 Jo<q<i/i (27r)2 

(26) 

Here we have neglected the off-diagonal elements of the 
Cooperon, as well as the g ^ corrections to the Hikami 
box in 1)25(1 . This is justified since those terms do not 
contain logarithms in x (or dephasing length) , in contrast 
to the diagonal Cooperon contributions in H2t)|l [see eqs. 
(EHIl,® below]. 



77 _ 



1 



^7 



(28) 



where for the left (L) and right (R) vertex parts we use 
the notations 



L{p) = Ga(p)pj; — Gr(p), 



(29) 



e 7 



R{-p) = Gr(-p) - -A,) Ga(-p) « ^g^^ 
~ Gr(-p) {-u^pf) Ga(-p). 

Retaining only leading order, we have approximated the 
charge current vertex in (|3l)|l in accordance with the dis- 
cussion at the end of Sec. IIIII [The same approximation 
is made in the calculation of the weak localization cor- 
rection to the charge conductivity^.] The quantities A 
and B in 1281) arc defined as 



= Tr{i(p)a^Gl(-p)a^}, (31) 

- Tr {<T^i?(-p) [GA(p)a^]'^} , (32) 

AJ^^ = Tr {L(p)a^ [a^GR(-p)]^} , (33) 

B^!'> = Tr {a^GR(p)a^ [R{-p)f} ■ (34) 
Due to the symmetry properties 



{Ai^r^^Ai'^, {Bfr^-Bf, 
(An*=^r, {Brr^Br, 7^2, 



(35) 



we see that for all 7 V^'^ = (V^^)* in 

In addition, for every diagram in Fig. [21 one needs to 
consider the corresponding renormalization of the spin 
current vertices, which turns out to be of the same order 
in {pfI)^^ as the bare L{p), 



L(p) =L(p)+^Ga(p)-7^Gr(p) 
xD^'''TT{a^'L{p)}. 



(36) 



The corresponding diagrams arc depicted in Fig.|31 Their 
contributions add to the ones of the bare diagrams in 
Fig. 12 Note that the symmetry properties 1)35(1 remain 
valid when the renormalization 1(36(1 is taken into account 
in the expressions for A and B. 

Similary, the charge current vertex (cx R(jp)) gets 
renormalized, but unlike before this renormalization is 
of higher order in {pfI)~^ [see at the end of Sec. II Vj . 
and hence can be neglected. Still, we note that like in 



5 




FIG. 3: The weak localization diagrams from Fig. I^jwith the 
renormalized spin current vertex (oc L{p)), see eq. l|36|l . The 
analogous renormalization of the current vertex (oc -R(p)) is 
of higher order in {pfI)~^ and can be neglected. 



Sec. IIVI one can demonstrate that this renormalization 
of i?(p) results in cancellation of the "anomalous" term 
(X Ax in the charge current vertex in H3U|I . 

Now we have to obtain an expression for J (2^)2 C'^'^ (q) 

in H26() . For this we can make use of the results of Refi2^, 
where an analogous approach has been used to study the 
weak localization correction to the charge conductivity in 
the same model system. The Cooperon derived in Refi^ 
is connected with C in eq. (|23() via 



n = in (|10() . Adding up all three terms of the Hikami 
box we obtain V^^{0) = in Thus, we see that 

the weak localization contribution to the spin Hall con- 
ductivity, i.e. CTyi^'' (see eq. (|26ll ). vanishes in the limit 
uj = 0. 

The evaluation of the zero- and one-loop diagrams per- 
formed in this and the preceding sections is valid up to 
the order l/ppl. If one would like to go beyond that or- 
der (which is not attempted here) , many more corrections 
(e.g. of the type l/E-pr) and many more diagrams had 
to be taken into account. For instance, one would also 
need to go beyond the self-consistent Born approximation 
for the self-energy and to retain crossing diagrams. Of 
course, the weak localization contribution which includes 
the Cooperon (resulting from time-reversed paths) is spe- 
cial in the sense that (apart form the log-dependence) 
it is sensitive to phase coherent effects probed by an 
Aharonov-Bohm flux. Thus, other contributions, even 
if they are of the same order in 1/ppl, could be ignored 
provided one would be interested in the phase sensitive 
contributions only. 



C^^ (q) 



1 



E 



where C^^ 



is the Cooperon from Refi^. Then, 
eqs. (13) and (14) in Refi^ we obtain that 



C'^''(q) 



ttPt 



:(/,/, log ^,/), 



(37) 



from 



(38) 



where** stands for the dephasing length, and 



VI. SPIN CURRENT AND SPIN PRECESSION 

We have seen that in perturbation theory the spin- 
Hall conductivity vanishes in the zero frequency limit, in 
leading (oc (1/pfO°) a?^rf subleading (cx (1/pfO^) order. 
This result suggests that the vanishing of the spin current 
is an exact property of the system under consideration. 

Indeed, we give now a simple argument to support this 
claim, based on the generalization of the Hamiltonian 



In-^i, for .T<T^, 



In- 



for 



< a; < 1. 



(39) 



We note that due to this logarithmic dependence on x the 
Cooperon itself can become large for sufficiently small x 
and large dephasing length L^. However, this has no con- 
sequences for the spin Hall conductivity eq. (|26|1 . since 
the Hikami box contributions add up to zero, as we shall 
see next. Indeed, after some lengthy but straightforward 
calculation we find for the diagonal terms of the Hikami 



^^'(P, r) = + a {diPy - a2Px) + ^^2) 

+/3 {criPx - cr2Py) + erE, 

where E is the applied static electric field, and (3 is the 
amplitude of Dresselhaus spin-orbit interaction. 

From the Heisenberg equation of motion for the spin 
s = a/2 of the electron, i.e. s(t) = ^s(t) = i[H'^s]{t), 
we obtain a relation between spin precession and spin 



curren 



.43 



yj'^ = (o,o,o,o)-HO(x^), 



impp 



2 ^3 



16 



-(x2,8- 



15x^ 



-x\8~ Ux^). 



(40) 



(41) 



Quite remarkably, the terms independent of x come from 
the vertex renormalization (2nd term on the rhs of eq. 
(|36|l '). where the only non-zero contribution arises from 
the diffuson matrix element D^^ oc \/x^. Each of the 
terms VJ, and Vc, when inserted into eq. H26|) . gives a 
contribution to (Jyx oc e, which corresponds to the term 



1 
2m 

_ 1 asx - (3sy 



Jy 



2m o? 
1 



/?2 



2m ^2 



(43) 



(44) 



The uniform spin current operator (z-component of the 
spin) is obtained from Q and given in first quantization 
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notation by = (l/2){v, s^}, with (spin-dcpcndcnt) 
velocity operator v = i[H' , r]. 

From the exact relations H44|l we see that the spin cur- 
rent can be experimentally accessed by measuring the 
spin precession or the change of magnetization in time, 
which can be done e.g. with optical methods^iS^iSk. 

In case of a = or /3 = H44|l provides a simple 
physical interpretation of the spin current: The spatial 
X (i/)-component of the spin current which carries the z- 
component of the spin is, up to a coupling constant, given 
by the total time derivative of the x (?;)-component of the 
spin. 

Next, let us consider the expectation value (including 
disorder average) of (|44|l for t — > cx). The presence of dis- 
order is expected to provide a relaxation mechanism such 
that the system can reach a steady nonequilibrium state 
(after some transients) when driven by an external elec- 
tric field (assumed to be constant in time). The weaker 
the disorder the longer it takes to reach this stationary 
limit, but for any finite amount of disorder (finite density 
of random impurities), it will be reached eventually. In 
particular, the magnetization (s)(i) of such a state also 
reaches a stationary value (possibly different from zero) 
and does not depend on time anymore. [Phenomeno- 
logically, the approach to such a state is described by a 
Blocli equation for the spin dynamics; see also below.] 
Thus, it immediately follows that the rate of magnetiza- 
tion change, (s)(t), must vanish in such a steady state, 
and, consequently, the spin current vanishes as well. In 
other words, under the stated conditions the uniform spin 
current vanishes in the long time limit if and only if the 
magnetization reaches a stationary state, i.e. 



lim(f=)(i) = 

t — ^oo 



lim (sfe)(t) = const. 



■ = x,y. 
(45) 

The result H45|l holds for any form of the static (spin- 
independent) impurity potential U{r) in the Hamiltonian 
(|42|l . including the special case of isotropic scattering po- 
tentials considered in previous sections. Also, the argu- 
ment is valid for any strength of the electric field, and, 
hence, applies to the special case of linear response con- 
sidered in the previous section [with {jy') = jy']- Thus, 
from = {jy'){t oo) = cTyxi^^ = 0)Ex we conclude 
that the spin Hall conductivity cry^{uj = 0) vanishes ex- 
actly (i.e. in all orders of the disorder potential) for any 
finite amount of disorder. 

In the special case^ a = /?, it is convenient to rewrite 
H' in (|42|l in the coordinate system, rotated counter- 
clockwise by 7r/4 in the xy-plane: 



i/R(p',r') = i/'(i?V4P,^V4r) 



^ 2aa'2P^ + U'{y' ) + er'E' = i^Ro+er'E' 

2m 



(46) 



where 



Rtt/4 



1 



(47) 



and for convenience we have introduced an alternative 
set of Pauli matrices 



^12 = ■^(0-2 ±(7l). 



0-3 • 



(48) 



Suppose before switching on the electric field (at t = 
0) the system was described by an equilibrium density 

matrix po = e~"^°^'^/Z. [It implies (j'^=)(t = 0) = 0.] 
Then the spin current in the rotated system is given by 



(j^-)W = Tr 



— — e ^ poe 
2 m 



(49) 



Since 



Tr [a'M = 0, Tr [a^.i/Ro] =0, k = 1,3, 

spin spin 



Tr 

spin 



2 m 



0, yt, 



(50) 



so that (j^-)(i) = for all t > 0. Thus, once (j'*-) is 
zero at t = 0, it remains zero forever due to the fact that 
(sy — Sx) is a conserved quantity for the Hamiltonian l|42|l . 

We consider now /3 = 0. The approach of (s)(t) to 
its stationary value (in the linear response regime) can 
be easily illustrated by taking the inverse Laplace trans- 
form of eqs. (jS)) and H21() (restricting ourselves to the 
Boltzmann value). This involves solving a cubic equa- 
tion for obtaining the poles with respect to uj. The re- 
sulting expression for {sy){t) consists of two parts (too 
lengthy to be written down here), one coming from the 
real pole and an oscillatory one coming from the pair of 
complex conjugate poles. For <IC 1 and i 3> r, only 
the first part is relevant, which has no oscillations and 
decays exponentially to zero: 



{§y){t) = -~2ma{j''){t) = -p-max'^E, 

in 



-t/T 



t > T, 

(51) 

where T/2 = (A^r) ^ is the well-known Dyakonov-Perel 
spin relaxation time^. For > 1 and t <t, the oscilla- 
tory part in {sy){t) becomes dominant, with period 1/A 
and exponential decay with rate 1/t; its time dependence 
for a particular value of x is illustrated in Fig. 01 

The above consideration can be generalized to the case 
with electron-electron interaction. In this case, the total 
spin of the system, '^^Si_k^ and the total spin current, 
X^iii*!' ^ ~ ^tU- obey the same equation as before, i.e. 
eq. H51|l . Thus, the same argument goes through as well, 
showing that the spin current must vanish in the sta- 
tionary limit also for interacting systems-provided this 
limit exists, which, again, we expect to be the case in the 
presence of any finite amount of disorder. 

Finally, in the absence of disorder, no stationary limit 
of {s){t) can be reached, i.e. the magnetization can 
change indefinitely, and thus there can be a finite spin 
current in this very particular case. [We note that 
the physically observable quantity, the magnetization 
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FIG. 4: The time evolution of {Sy{t)) for a; = 4. The period of 
oscillation is 1/A, and the exponential decay time is r. Note 
that for any x, {sy{t)) = at t = 0. 

change, vanishes for vanishing Rashba coupling.] How- 
ever, when the spin current approaches a non-zero but 
constant value for t ^ oo, as obtained in the linear re- 
sponse regime for a clean systemic, the magnetization 
(sfc)(t) actually grows linearly in time in the asymptotic 
regime t ^ oo. This, of course, shows a breakdown of 
the linear response approximation in this case since the 
bound for a spin 1/2, i.e. \{sk){t)\ < 1/2, is violated. The 
terms beyond linear response would restore the bounded 
and oscillatory behavior of {sk){t). 



VII. CONCLUSIONS 

In conclusion, we have calculated the Boltzmann and 
weak localization contributions to the spin-Hall conduc- 
tivity (7y^ in an infinite-size 2D disordered system with 
Rashba spin-orbit interaction. In our calculation we 
had to consider the weak localization contribution to the 
, together with the Boltzmann contribution, since dia- 
grams belonging to these two contributions have the same 
order of magnitude. In the diffusive and static limit, the 
spin-Hall conductivity vanishes. We have shown that the 
spin current is given by the time derivative of the mag- 
netization change and have given general arguments why 
the spin current must vanish for any finite amount of 
disorder that permits to reach a stationary limit of the 
system when driven by a constant electric field. It will be 
interesting to see how far our arguments can be extended 
to other types of spin orbit interaction. 
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